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摘要 : ADER Oedaleus infernalis Saussure 是 一 种 在 我 国 分 布 广泛 、 对 农 牧 业 生产 危害 严重 的 经 济 害虫 。 本 文 
XR NZEBR 10 个 地 理 种 群 的 线粒体 COI 基因 序列 进行 测序 和 分 析 , 利用 DnaSP 5.0 和 Arlequin 3.5. 1.2 软件 对 
该 蝗虫 种 群 间 的 遗传 多 样 性 、 遗 传 分 化 程度 、 基 因 流 水 平 及 分 子 变 异 进行 了 分 析 , 建立 了 单 倍 型 贝 叶 斯 系统 发 育 
进化 树 和 单 倍 型 网 络 图 。 结 果 表 明 : 在 所 分 析 的 144 个 序列 样本 中 , 共 检 测 到 21 种 单 倍 型 , 其 中 1 种 单 倍 型 为 10 
个 地 理 种 群 所 共享 。 总 群体 的 单 倍 型 多 样 性 指数 为 0.653, 各 地 理 种 群 单 倍 型 多 样 度 范围 在 0.423 ~ 0. 790 之 间 。 
总 群体 和 各 种 群 的 Tajima’ s D 检验 结果 皆 不 显著 ,说 明 该 种 害虫 在 较 近 的 历史 上 未 经 历 群 体 扩 张 。 总 群体 的 遗 
传 分 化 系数 Gst X 0.04436, 固定 系数 Fst 为 0.05255, 基因 流 Nm 为 9.01。AMOVA 分 子 方差 分 析 结 果 表 明 , RIA 
小 车 蝗 的 遗传 分 化 主要 来 自 种 群 内 部 , 种 群 间 的 遗传 变异 水 平 较 低 。 各 地 理 种 群 的 遗传 距离 的 大 小 与 其 地 理 距离 
间 没 有 显著 的 相关 性 。 贝 叶 斯 系统 发 育 进 化 树 与 单 倍 型 网 络 图 显示 , 黄 用 小 车 星 各 地 理 种 群 中 的 单 倍 型 散布 在 不 
同 的 分 布 群 中 , 分 布 格局 较为 混杂 , 未 形成 明显 的 系统 地 理 结构 。 研 究 结 果 揭 示 , 黄 肥 小 车 星 各 种 群 间 的 基因 交 
流 并 未 受到 地 理 距 离 的 影响 。 
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Analysis of genetic differentiation and gene flow among different 
geographic populations of Oedaleus infernalis ( Orthoptera: Acrididae ) 


based on mtDNA COI gene sequences 

SUN Wei, ZHANG Zhu-Ting, DONG Hui, QIAN Hai-Tao, SHI Yu, XIE Li-Na, CONG Bin ( College 
of Plant Protection, Shenyang Agricultural University, Shenyang 110161, China) 

Abstract: Oedaleus infernalis Saussure is one of the most important agricultural and grass forage pests, 
which is widely spread in China. The partial sequences of mtDNA COI gene in 10 geographic populations 
of O. infernalis were sequenced and analyzed. The genetic diversity, genetic differentiation, gene flow 
and molecular variance were analyzed using DnaSP 5. 0 and Arlequin 3. 5. 1. 2. Bayesian inference ( BI) 
phylogenetic tree and the haplotype network of COI haplotypes were constructed. The results showed that 
a total of 21 haplotypes were defined among 144 gene sequences, including one haplotype shared by ten 
populations. The total haplotype diversity of all populations was 0. 653 and the range of haplotype 
diversity of different populations was 0. 423 —0. 790. The Tajima’ s D test indicated that there might be 
no population expansion for the total in recent history. The total Gst was 0. 04436, Fst was 0. 05255 and 
the gene flow Nm was 9. 01 among 10 populations. Molecular variance analysis ( 4AMOVA) demonstrated 
that observed genetic differentiation mainly occurs within populations and just a little differentiation exists 
among populations. The genetic distance among populations was not correlated with the geographic 
distribution. The Bayesian inference ( BI) phylogenetic tree and haplotype network showed that 
haplotypes were distributed in different clades and no obvious geographical structure formed. The results 
suggest that the gene flow among O. infernalis populations is not affected by the geographical distance. 
Key words: Oedaleus infernalis; geographic population; mtDNA COI gene; haplotype; genetic 


differentiation; gene flow 
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3a He/] 4E Oedaleus infernalis Saussure 在 系统 
TXR ERJ T EW OH (Orthoptera), 8 P} 
( Acridoidea ) , "i EL ( Acrididae ), PE XH 8E Wy $ 
( Oedipodinae) , /| 4:1 gj Oedaleus, 在 我 国 主要 分 
fi PF. HM. UUT, TR, HON. EAE 1l 
西 、 河 北 、 山 东 等 省 区 ， 国 外 的 前 苏联 、 日 本 、 驼 
十、 韩国 也 有 分 布 。 该 虫 在 北京 以 北 地 区 每 年 发 生 
1 代 , 华北 中 、 南 部 和 中 原 地 区 2 代 , 均 以 卵 块 在 
十 中 越冬 。 黄 肥 小 车 昌 属 杂食 性 害虫 , 喜 食 插 草 、 
针 碟 、 冰 草 、 小 麦 、 玉 米 、 高 梁 等 开本 科 牧 草 及 作 
物 , 还 可 取 食 大 豆 、 马 铃 昔 等 作物 , 是 我 国 重要 的 
ABO E RS Reo NAE RE SERES] iz 
分 布 ( 李 宏 实 和 张 喜 军 ，1991; 任 炳 忠 , 2002; FW 
5%, 2007), 是 我 国 北方 草原 的 主要 优势 种 蝗虫 ， 
对 畜牧 业 的 可 持续 发 展 造成 严重 影响 ; 在 河北 、 山 
东 等 省 部 分 农 区 黄 肥 小 车 蝗 的 发 生 面 积 大 , 为害 程 
度 重 , 是 威胁 小 麦 及 大 秋 作 物 生产 的 主要 害 忠 之 一 
( 任 春光 , 2009; 田 方 文 和 李 振 国 , 2009), 

HAI FEXR TEINER EYRE 
分 布 、 危 害 损失 估计 和 分 子 系统 学 等 方面 进行 了 广 
泛 深入 的 研究 ( 何 建 平 等 ,2006; Xims, 2008; 
Fartmann et al., 2008; 齐 巧 丽 等 ,2009 ; Yoshioka et 
al., 2010; Kindler et al., 2012) , 但 迄今 有 关 黄 腑 小 
车 临 地 理 种 群 的 遗传 结构 及 分 化 的 报道 较 少 。 相 对 
FKE, WWE tE Oedaleus decorus asiaticus Bei 
Bienko 等 报道 有 迁 飞 习性 的 蝗虫 , WEIR 
为 是 多 进行 近 距 离 扩 散 克 食 、 求 偶 交 配 或 逃逸 ,不 
远 距离 迁 飞 的 土 螺 。 由 于 食性 杂 、 分布 广 , 不 同 地 
理 种 群 的 黄 肥 小 车 旦 在 生物 学 特性 及 生活 史上 有 所 
差异 。 环 境 因子 的 不 对 称 性 ， 是 否 造 成 了 基因 流 的 
不 随机 性 ( Avise,1992; FE, 2010), fr HE 
车 星 地 理 种 群 间 出 现 不 同 程度 的 遗传 变异 ,从 而 构 
成 星 害 日 益 严 重 的 内 在 遗传 因素 ,这 一 科学 问题 有 
竺 解决。 

随 着 分 子 生物 学 的 发 展 和 应 用 , 分 子 标记 技术 
越 来 越 多 地 被 应 用 到 蚂 虫 种 群 遗 传 研究 中 去 ， 如 
Cooper 等 (1995 ) 用 核 DNA 序列 标记 研究 了 欧洲 临 
Hi Chorthippus parallelus ( Zetterstedt ) 冰期 后 重建 相 
天 的 遗传 亚 种 群 模式 ; Zhang SE (2009 ) 用 微 卫 星 分 
子 标记 对 中 国 发 晶 种 下 阶 元 进行 了 重新 划分 ,并 探 
讨 了 飞 蝗 的 历史 演化 过 程 。 在 众多 的 分 子 标记 技术 
中 , 线粒体 DNA. (mtDNA ) 因 其 操作 简单 、 母 系 遗 
传 、 进 化 速率 快 、 找 贝 数 多 等 特点 ( Moritz et al., 
1987; Harrison, 1989; Wolstenholme，1992 ) ,被 广 


泛 应 用 于 蝗虫 的 系统 进化 及 种 群 的 遗传 变异 的 研究 
中 ( 蒋 国 芳和 郑 哲 民 , 2001; 叶 海燕 和 黄 原 , 2011) 。 
本 实验 选用 线粒体 COI 作为 分 子 标记 , 对 黄 胜 小 车 
蝗 10 个 地 理 种 群 的 遗传 结构 和 基因 交流 模式 进行 
研究 ,分 析 其 不 同 地 理 种 群 的 遗传 分 化 程度 及 遗传 
多 样 性 , 明确 其 生态 适应 性 ,探讨 种 群 动态 及 星 灾 
发 生 的 内 在 遗传 因素 ,为 针对 该 蝗虫 有 效 防 治 策略 
的 制定 提供 分 子 生物 学 的 基础 理论 资料 。 


1 材料 与 万 法 


1.1 供 试 虫 源 

本 研究 中 选用 的 10 AP 38 Ee 7] 4 EH REB RF, 
T 2011 年 采 自 吉林 、 辽 宁 、 内 蒙古 、 河北 、 陕 西 、 
山西 、 宁 夏 7 个 省 份 。 各 种 群 的 具体 采集 信息 及 地 
理 群 体 缩写 见 表 1。 将 采集 到 的 成 虫 浸泡 于 无 水 乙 
醇 中 , 浸泡 的 第 一 个 月 需 更 换 无 水 乙醇 3 ~4 次 , 后 
放置 于 -70Y 超 低温 冰箱 中 长 期 保存 待 用 ( 张 民 照 
和 康乐 , 2001; 张 建 珍 等 , 2004 ) 。 
1.2 基因 组 总 DNA 的 提取 

基因 组 DNA. 的 提取 参照 刘波 等 (1999 ) RB I 
华 等 (2011 ) 报道 的 方法 并 稍 加 改进 。 取 适量 晶 虫 
后 足 股 节 肌 和 肉 ， 冲 洗 后 剪 碎 。 加 入 400 pL 196 
CTAB 533 , i je 30 s, 377C 恒温 水 浴 1 ho MEH 
酶 K 0.06 mg 和 10% SDS 40 uL, 65*C 恒温 水浴 处 
JE 4 ho WMA 2 倍 体积 氯仿 : ERE (24: 1, v/v) , 充 
分 混 匀 ,10 000 r /min 离心 10 min。 取 上 清 液 加 等 
体积 异 丙 醇 ，- 20 沉淀 DNA 2 h。70% 乙醇 洗涤 
后 干燥 。 加 TE 缓冲 液 及 适量 RNase 37% 恒温 水 浴 
保温 1 h。 琼 脂 糖 凝 胶 电泳 及 Nanodrop 2000c 4X di 
( Thermo Scientific 公司 ) 检测 DNA 浓度 和 纯度 , S 
格 后 -20Y 保存 备用 。 
1.3 了 PCR 扩 增 和 产物 测序 

扩 增 的 目的 片段 为 Lunt 等 (1996 ) 15:38 E] Eb R 
线粒体 DNA COI 中 一 段 进化 速率 较 快 、 可 用 于 种 
下 遗传 分 化 研究 的 基因 序列 。 上 游 引 物 为 UEA3: 
5'-TATAGCATTCCCACGAATAAATAA-3', 下 游 引 物 
Jj UEA4: 5'-AATTTCGGTCAGTTAATAATATAG-3' 
( Lunt et al., 1996; Zhang and Hewitt, 1997) , 

PCR 扩 增 反应 体系 为 50 pL, 其 中 模板 2 pL 
( 约 50 ng) 、 上 下 游 引物 (10 pmol/L) && 0. 75 pL, 5 
uL 10 x PCR buffer, 3.8 uL MgCl, (20 umol/L), 1 
uL dNTPs (10 iumol/L) , 1 uL Taq 酶 (2 U/uL), 
超 纯 水 补足 。 扩 增 条 件 为 94Y HRAETE 4 min; 94'C 
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表 1 黄 肥 小 车 蝗 采 集 信息 
Table 1 Collecting data of Oedaleus infernalis in this study 


种 群 代 码 采集 地 点 

Population code Collecting locality 
FK 辽宁 法 库 Faku, Liaoning 
BP 辽宁 北 票 Beipiao，Liaoning 
NA 吉林 农安 Nong an, Jilin 
TL 内 蒙古 通辽 Tongliao, Inner Mongolia 
CR 内 蒙古 格 日 朝 鲁 Gerichaolu, Inner Mongolia 
XL 内 蒙古 锡林浩特 Xilinhaote, Inner Mongolia 
CZ 河北 沧州 Cangzhou，Heibei 
HL THW Helan, Ningxia 
QX 山西 清 徐 Qingxu, Shanxi 
WQ 陕西 吴起 Wuqi, Shaanxi 


变性 40 s, 43% 退火 1 min, 72% 延伸 1 min 20 s, 4 
计 35 个 循环 ; 最 后 727C 延伸 10 mins PCR 产物 4 
pL T 1.296 的 琼脂 糖 凝 胶 进 行 电 泳 检测 , 确定 扩 增 
成 功 的 委托 北京 易 国 生物 技术 有 限 公 司 测 序 。 
1.4 数据 统计 与 分 析 

测序 结果 利用 Chromas 软件 读 取 、 观 察 峰 信 并 
用 DNAStar 软件 包 对 序列 进行 编辑 校正 。 将 处 理 后 
的 序列 在 NCBI 网 站 上 进行 BLAST 比 对 、 确认 。 通 
过 DNAMAN 软件 进行 多 重 序 列 比 对 。 使 用 
MEGA4.0 软件 (Tamura et al., 2007 ) 分析 碱 基 组 
成 、 变 异 位 点 、 简 约 信息 位 点 及 碱 其 转换 丰 换 比 
值 ; 基于 Kimura-2-Parameter 模型 计算 单 倍 型 之 间 
的 遗传 距离 ; 以 网 翅 蝗 科 绿 牧草 蝗 Omocestus 
viridulus ( L.) COIL 基因 片段 作为 外 群 , 采 用 
MrBayes v3. 1. 2 软件 (Huelsenbeck and Ronquist, 
2001) 构建 贝 叶 斯 ( Bayesian inference, BI) |, $& 7) 
文 置信 和 度 以 后 验 概 率 来 评价 。 单 倍 型 网 络 图 由 
NETWORKA. 6. 1. 0 软件 (Bandelt et al., 1999) 
绘制 。 

DnaSP 5.0 ( Librado and Rozas, 2009) 软件 分 析 
10 个 种 群 内 的 单 倍 型 多 样 度 Hd, EREE Pi, 
核 甘酸 平均 差异 数 K, 种 群 间 核 车 酸 平 均 差 异 数 
Key, CH BRI SCIRE Dxy、 遗 传 分 化 系数 Gst 与 固定 
系数 Fst 等 分 子 遗 传 学 系数 ， 并 进行 中 性 检验 
(Tajima，1989 ) 。 应 用 Arlequin 3. 5. 1. 2 软件 
( Excoffier et al., 2005) 进行 AMOVA 分 子 变异 分 析 ， 
计算 种 群 间 的 遗传 距离 。 依 据 经 纬度 计算 各 采集 地 
之 间 的 地 理 距离 。 使 用 GenAIEx 6. 41 软件 (Peakall 
and Smouse, 2006) 检验 种 群 间 的 遗传 距离 与 地 理 





地 理 坐 标 样本 量 采集 时 间 
Geo-coordinates Sample size Sampling date 

42.54?N, 123.31?E 15 2011.8 
41.91? N, 120. 79?E 13 2011.10 
44. 34? N, 124. 82?E 15 2011.8 
43. 60?N, 122.29?E 15 2011.8 
44. 49?N, 120. 05?E 15 2011.8 
44.21?N, 115. 98?E 13 2011.9 
38.29?N, 116. 89°E 15 2011.9 
38. 74?N, 106. 10? E 14 2011.9 
37.67?N, 112.25?E 15 2011.9 
37.07?N, 108. 18? E 14 2011.9 
距离 的 相关 性 。 
2 ”结果 与 分 析 


2.1 碱 基 组 成 及 序列 变 

对 10 个 地 理 种 群 共 144 3E AP BRL RS ^] ita 
COI 基因 片段 测序 , 经 多 重 序列 比 对 , 截 去 测序 误 
差 较 大 的 区 域 , 得 到 长 度 为 305 bp 的 片段 ,与 
GenBank 上 已 发 表 的 黄 肥 小 车 曙 COI 序列 
( GenBank 登录 号 AB497592. 1 ) 一 致 性 最 高 可 达 
99% , 确认 为 目的 片段 。 在 COI 基因 片段 序列 中 ， 
没有 发 现 碱 基 的 插入 和 缺失 , 保守 位 点 290 个 。 多 
态 性 位 点 15 个 ， 占 所 测序 列 长 度 的 4.92% , 其 中 
简约 信息 位 点 8 个 , 自 诊 位 点 7 个。 检测 到 的 总 突 
变数 为 16, 其 中 发 生 转 换 的 位 点 数 11 个 , REN 
换 的 位 点 数 5 个 ,转换 以 CT HE, BREL CA 为 
E, 1 个 ACT 三 碱 基 变 异 位 点 ， 转换/ 颠 换 侦 倚 
率 R 为 2.9。 从 已 报道 的 研究 中 可 知 ， 杀 缘 关 系 较 
近 的 分 类 阶 元 中 核 音 酸 蔡 换 的 转换 多 于 斯 换 ， 而 远 
缘 的 分 类 阶 元 中 斯 换 多 于 转换 , 表现 出 明显 的 差异 
(Simon et al., 1994; 李靖 等, 2010) , 本 研究 结果 与 
其 相近 。 在 COI 基因 片段 的 组 成 中 , 4 种 碱 基 的 平 
均 含 量 依 次 为 A: 36.1996, T. 31. 7296, C: 
18.9896, G: 13. 09976, A + 了 T 的 平均 含量 为 
67.91% , 表现 出 A/T Bà 3E 2H X fi I] TE ,— — C. ATL ER 
虫 mtDNA 序列 A +T 碱 基 舍 量 较 高 的 报道 相符 合 
(Liu and Beckenbach, 1992) , 
2.2. 单 倍 型 分 析 

实验 中 获得 的 黄 肥 小 车 蝗 10 个 地 理 种 群 的 21 
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种 单 倍 型 ， 记 录 为 HI ~ H21 (GenBank 登录 号 : 
KC297197 ~ KC297217) 。 应 用 Kimura-2-Parameter 
模型 计算 单 倍 型 之 间 的 遗传 距离 ,结果 显示 , 不 同 
单 倍 型 遗传 距离 在 0. 003 ~ 0. 027 之 间 , 平均 遗传 
EZ J 0.01019, 

TE Ir x PLE RES] CO 基因 序列 的 21 种 
单 倍 型 中 , 单 倍 型 Hl 出 现 频 率 最 高 (83/144 ) ， 为 


H20 


共享 单 倍 型 ， 占 所 有 检测 个 体 的 57.6490; 其 次 出 
现 频率 较 高 的 单 倍 型 是 H2， 占 总 个 体 的 9. 7290 
(14/144) , 出 现在 所 检测 的 6 个 种 群 中 ; HIS 单 倍 
型 占 样 本 的 6. 94% (10/144) ， 全 部 出 现在 农安 
(NA)、 法 库 (fK) 两 个 地 理 种 群 中 。 由 单 倍 型 网 络 
图 (图 1) 可 清晰 地 看 出 各 单 倍 型 在 各 地 理 种 群 的 分 
布 情况 及 单 倍 型 之 间 的 演化 关系 。 
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图 1 RANE mtDNA COI 的 基因 单 倍 型 网 络 图 
Fig. 1 Median-joining haplotype network of Oedaleus infernalis based on COI gene of mtDNA 
圆 面积 代表 单 倍 型 出 现 频率 , 彩色 肩 形 面积 代表 各 样品 种 群 在 同一 单 倍 型 中 所 占 的 比例 。Cirele areas are proportional to haplotype frequencies, 


while colored portions represents the proportions of the same haplotype that occurs in each sampling region. 


2.3 单 倍 型 系统 进化 分 析 

图 2 为 采用 贝 叶 斯 法 (4BI) AEKA 
21 种 单 倍 型 与 外 群 之 间 的 系统 发 生 树 。 结 果 显 示 ， 
ETRE COI 基因 片段 的 各 单 倍 型 均 明显 
与 外 群 分 开 , 但 单 倍 型 之 间 并 未 显示 出 明显 的 地 理 
谱系 ， 各 分支 置信 和 度 也 较 低 (后 验 概率 低 于 50% 
HJ, 系统 树 未 作 显 示 ) 。 尽 管 有 距离 较 近 的 地 理 种 
群 的 单 倍 型 优先 聚 在 一 个 分 文 (如 GR. 的 H5 A TL 
种 群 的 H9) , 但 也 有 地 理 距 离 相 距 较 远 的 单 倍 型 聚 
在 一 起 的 情况 (如 CR 的 HI2 和 QX 的 HI4) 。 总 体 
来 看 , 单 倍 型 的 拓扑 结构 并 未 反映 出 与 地 理 位 置 相 
关 的 信息 , 未 形成 明显 的 系统 地 理 结 构 。 


2.4 遗传 多 样 性 

由 表 2 所 示 , 10 个 地 理 种 群 的 总 群体 单 倍 型 多 
样 度 Hd 为 0.653, 核 苦 酸 多 样 度 Pi H 0.00364 , 核 
苷 酸 平均 差异 数 开 为 1.10975。 各 地 理 种 群 单 倍 型 
数量 范围 在 4 ~8 之 间 , 平均 具有 单 倍 型 4.9 个 , 其 
中 格 日 朝 鲁 (CR ) 种 群 单 倍 型 数量 最 丰富 , 在 15 个 
检验 个 体 中 共和 定义 了 8 种 单 倍 型 。 不 同 地 理 种 群 
COI 基因 的 单 倍 型 多 样 度 Hd 范围 为 0. 423 ~ 
0.790, 平均 值 为 0.612, 具有 最 多 单 倍 型 数 的 格 日 
朝 鲁 ( GR) 种 群 内 单 倍 型 多 样 度 最 高 , 北 票 (BP)、 
锡林浩特 (XL ) 种 群 内 的 单 倍 型 多 样 度 最 低 。 各 种 
群 平均 核 苷 酸 多 样 度 Pi J 0.00341, 不 同 种 群 的 Pi 


$^ 


Omocestus viridulus 





图 2 MIAA CO 单 倍 型 系统 树 
Fig. 2 Bayesian inference ( BI) phylogenetic tree of COI 
haplotypes of different Oedaleus infernalis populations 
VAS Bo AE HOA ARRE, Ar sc E BORSE OS Jo T LOS fo Omocestus 
viridulus was used as the outgroup. Numbers on the tree represent the 


values of posterior probability values. 
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变化 范围 在 0. 00179 ~ 0.00556 之 间 。 核 车 酸 平均 
差异 数 天 在 各 地 理 种 群 的 分 布 范围 为 0. 54545 ~ 
1.69524, 均值 为 1. 03932, Pi. k 值 最 高 及 最 低 的 
种 群 分 别 是 格 日 朝 鲁 (GR) 和 次 兰 (HL) 种 群 。 各 地 
理 种 群 Tajima' s D 中 性 检测 均 无 显著 性 差异 , ERAS 
兰 (HL)、 锡 林 (XL) 两 个 种 群 为 0. 10 > P >0.05, 
其 他 种 群 丝 为 P >0.10。 总 群体 的 Tajima' s D 检测 
为 负 值 , 不 显著 (0.10 >P >0.05)。 
2.5 遗传 分 化 分 析 

应 用 DnaSP 5.0 软件 计算 反映 黄 肥 小 车 星 种 群 
间 基 因 差 异 程度 的 遗传 学 参数 ( 表 3, 表 4)。 结 
显示 , 各 地 理 种 群 总 遗传 分 化 系数 Gst 为 0. 04436 ， 
总 固定 系数 Fst 为 0.05255 ， 总 基因 流 Nm 为 9.01。 
各 种 群 间 核 甘酸 差异 数 Kxy 在 0. 5664 ~1.5956 之 
E, HEX 1.0967; 2P RE XE Dry Æ 0.0019 ~ 
0.0052 之 间 , 均值 为 0.0036; 遗传 分 化 系数 Gsi 在 
-0.0260 ~0. 1147 之 间 , 均值 为 0.0224; 固定 系数 
Fst 在 -0.0537 ~0.1804 之 间 , 均值 为 0. 0479 , 除 
农安 CNA) 种 群 与 部 分 种 群 遗 传 分 化 程度 达到 中 等 
水 平 (Weight,1978 ) ， 其 余 均 较 低 。 


HANERE [e] 3838 RR RE COI 单 倍 型 、 核 彰 酸 多 样 性 分 析 及 Tajima' s D 中 性 检验 


Table 2 COI haplotype diversity, nucleotide diversity and Tajima’ s D in different 


geographic populations of Oedaleus infernalis 





FCR IECDEBI2: 
倍 型 总 数 倍 型 多 样 度 AREFE 
种 群 年 作弄 PER l ^ ^ MN S Average number of 中 性 检测 中 性 检验 显著 性 
Number of haplotypes Haplotype diversity Nucleotide diversity 
Population , nucleotide differences Tajima’s D Statistical significance 
(H) ( Hd) (Pi) 
(K) 

BP 4 0. 423 0.00437 1. 33333 — 1.55085 P »0.10 

CZ 5 0. 733 0. 00306 0. 93333 —0. 77356 P »0.10 

NA 5 0. 752 0.00375 1. 14286 -0.22755 P »0.10 

GR 8 0. 790 0.00556 1. 69524 — 0. 76543 P »0.10 

HL 4 0.491 0.00179 0. 54545 — 1.59996 0.10» P »0.05 
FK 6 0. 758 0.00351 1. 07190 — 0. 23482 P »0.10 

QX 5 0. 562 0.00350 1. 06667 — 1.03336 P »0.10 

TL 4 0. 543 0.00281 0. 85714 -0.97211 P »0.10 

XL 4 0. 423 0.00202 0.61538 — 1.77497 0.10» P »0.05 
WQ 4 0. 648 0.00371 1.13187 — 0. 32983 P »0.10 
dd 21 0. 653 0.00364 1. 10975 — 1.64591 0.10» P »0.05 


Total population 
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Table 3 Pairwise average number of nucleotide difference Kxy (above the diagonal) and nucleotide divergence Dxy 


BP 
BP 

CZ 0.0039 
NA 0. 0045 
CR 0.0052 
HL 0.0030 
FK 0. 0043 
QX 0. 0040 
TL 0. 0037 
XL 0. 0032 
WQ 0.0038 

Table 4 
BP 

BP 

CZ 0.0349 
NA 0.1058 
CR 0.0505 
HL -0.0176 
FK 0.0743 
QX 0. 0043 
TL 0. 0248 
XL -0.0100 
WQ - 0.0537 


CZ 
1.1744 


. 0042 
. 0045 
. 0026 
. 0034 
. 0032 
. 0029 
. 0026 


O o O O oO o o Qc 


. 0036 


oococococc 


NA 


. 3846 
. 2667 


0052 
0031 
0038 
0042 
0039 
0031 


. 0043 


GR 
1. 5949 
1.3689 
1. 5956 


0038 
0050 
0046 
0042 
0040 


oooccococso 


. 0050 


( below the diagonal) among different populations 


HL FK QX TL XL WQ 
0. 9231 1.2992 1. 2051 1.1231 0. 9645 1. 1703 
0. 8061 1.0518 0. 9822 0. 8756 0. 7795 1. 0905 
0. 9576 1.1630 1. 2800 1. 2000 0.9590 1.3191 
1. 1546 1.5259 1.3911 1.2667 1.2205 1. 5381 

0.9293 0. 8303 0.7333 0. 5664 0. 8766 
0. 0031 1. 1037 1.0222 0. 8803 1. 2064 
0. 0027 0.0036 0. 9333 0. 8461 1. 1095 
0. 0024 0.0034 0.0031 0.7436 1. 0524 
0. 0019 0. 0029 0. 0028 0. 0032 0. 8956 
0. 0029 0. 0040 0. 0036 0. 0034 0. 0029 


表 4 种 群 间 遗 传 分 化 系数 Gst( 上 三 角 ) 与 固定 系数 Fst( 下 三 角 ) 


CZ 
0. 0466 


0.1804 
0.0399 
0.0827 
0.0468 
-0.0181 
-0. 0225 
0. 0066 
0. 0531 


O 0 O O oO oO qo 


NA 


. 1147 
. 0640 


. 1106 
. 1184 
. 0478 
. 1369 
. 1667 
. 0833 
. 1377 


GR 
0. 0269 
0. 0017 
0. 0372 


0. 0322 
0. 0933 
0. 0073 
— 0. 0075 
0. 0534 
0. 0810 


Pairwise Gst ( above the diagonal) and Fst ( below the diaognal) values among different populations 


HL FK QX TL XL WQ 
-0.0180 0.0513 - 0. 0052 0. 0045 - 0. 0178 0. 0026 
0. 0370 0.0018 0.0178 - 0. 0055 0.0427 — 0. 0002 
0. 0830 0. 0004 0. 0812 0. 0913 0.1079 0. 0721 
0.0110 0.0048 0. 0028 0. 0066 0.0269 0. 0031 

0.0347 — 0. 0152 - 0. 0007 — 0. 0260 0. 0130 
0.1298 0.0238 0.0149 0.0452 0.0135 
0. 0292 0.0312 - 0. 0093 - 0. 0052 0. 0037 
0. 0437 0.0565 - 0. 0306 0. 0045 - 0. 0043 
- 0. 0247 0.0417 0. 0061 0. 0098 0. 0224 
0. 0433 0. 0866 0. 0092 0. 0550 0. 0245 


AMOVA 分 子 变 异 分 析 结 果 显 示 ( 表 5)，, XR 
小 车 曙 10 个 地 理 种 群 间 的 遗传 变异 (5.33% ) i 远 小 
于 种 群 内 的 遗传 变异 (94.67% )。 膏 明黄 肥 小 车 曙 
的 遗传 变异 主要 来 自 种群 内 部 , 种 群 间 的 遗传 变异 
较 低 。 

依据 Kimura-2-Parameter 模型 基 


于 COI 3E pA FE 


(Slatkin, 1995) 。 由 表 6 可 知 , 各 种 群 间 遗传 距离 
范围 介 于 0 ~ 0. 2206 之 间 。 应 用 Mantel 相关 性 检 
验 检测 种 群 之 间 的 遗传 分 化 程度 是 否 与 地 理 隔 离 有 
关 。 结 果 表 明 , 成 对 种 群 的 遗传 距离 矩阵 与 采集 地 
点 地 理 距离 的 自然 对 数 和 矩阵 之 间 未 呈现 出 显著 相关 
性 (R=0.058, P 20.360 >0.05), 地 理 隔离 并 不 能 


列 计算 黄 肥 小 车 旦 不 同 地 理 种 群 之 间 的 遗传 距离 。 解释 种 群 间 遗 传 距离 的 差别 。 
RS EHE NES 10 个 种 群 的 分 子 变异 分 析 
Table 5 Analysis of molecular variance (AMOVA) among 10 populations of Oedaleus infernalis 
变异 来 源 自由 度 平方 和 方差 组 分 方差 比率 (% ) 
Source of variation df Sum of squares Variance components Percentage of variation 

种 群 间 Among populations 9 8. 660 0. 02995 Va 5.33 
种 群 内 Within populations 134 71.234 0.53160 Vb 94. 67 

总 变异 Total 143 79. 894 0. 56155 100 


8 期 孙 
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ARE] ENS Ep EPR 3038 ER A ( km) 自然 对 数值 (上 三 角 ) 与 遗传 距离 (下 三 角 ) 


Table 6 The natural logarithm of geographical distance ( km) (above the diagonal) and pairwise genetic distance 


( below the diagonal) among Oedaleus infernalis populations 


BP CZ NA GR HL FK QX TL XL WQ 
BP 6.26 6.05 5.68 7.18 5.43 6.74 5.42 6.15 7.04 
CZ 0. 0387 6.85 6.61 6.85 6.59 5.97 6.61 6.50 6.53 
NA 0. 1202 0. 2206 5.94 7.44 5.44 7.14 5.39 6.56 7.35 
GR 0.0519 0.0416 0. 1248 7.21 5.85 6.88 5.32 5.79 7.13 
HL 0 0.0801 0. 1181 0.0184 7.34 6.33 7.30 6.96 5.87 
FK 0. 0843 0. 0483 0.0511 0. 1083 0.1269 6.97 5.00 6.44 7.22 
QX 0. 0054 0 0.1593 0.0075 0.0212 0.0323 6.95 6.64 5.70 
TL 0. 0283 0 0. 2002 0 0. 0388 0.0583 0 6.24 7.20 
XL 0 0. 0040 0.0850 0. 0494 0 0.0355 0. 0033 0. 0081 6.92 
WQ 0 0.0574 0.1601 0.0871 0.0365 0. 0964 0.0097 0.0598 0.0232 
et al., 2012) 。 本 研究 结果 显示 , EUNT RE ADR 
3 讨论 理 种 群 内 具有 较为 丰富 的 单 倍 型 多 样 性 , VES] SCA 


细胞 色素 C. 氧化 酶 亚 基 ICOI) 基 因 是 13 个 重 
日 质 编码 基因 之 一 , 其 产物 细胞 色素 氧化 酶 亚 基 T 
是 细胞 色素 氧化 酶 的 重要 组 成 部 分 。 前 人 人 研究 表 
H, 该 基因 厂 段 不 同位 置 的 进化 速率 有 所 不 同 , 可 
通过 选取 适宜 的 基因 片段 , 进行 昆虫 高 级 阶 元 系统 
发 育 或 者 近 缘 种 、 种 下 分 类 单元 的 系统 关系 研究 
( Lunt et al., 1996; Zhang and Hewitt, 1997) 。 本 研 
究 选 取 的 是 CO 基因 序列 中 进化 速率 较 快 、 在 昆虫 
种 下 遗传 应 用 较 多 的 一 段 基 因 序 列 (Funk et al., 
1995; Bensasson et al., 2000) , XI me] AE k ST Hl 
理 种 群 遗 传 分 化 进行 了 分 析 。 

研究 绪 末 表明 ，10 NEJEDE a e E 
生 的 21 种 单 倍 型 中 , 单 倍 型 Hl 是 各 种 群 共 盏 的 主 
体 单 倍 型 , 可 认为 该 单 倍 型 是 较为 原始 的 、 能 够 适 
应 环境 变化 并 在 种 群 中 稳定 存在 的 优势 单 倍 型 ( 张 
立 勋 等 , 2005; 李 继 变 和 任 竹 梅 , 2009 ) 。 而 独立 存 
在 于 各 地 理 种 群 的 独 享 单 倍 型 则 说 明 各 地 理 种 群 在 
存在 着 一 定 基因 交流 的 同时 , 也 具有 一 定 程度 上 的 
遗传 分 化 。 综 合 单 倍 型 系统 树 及 中 介 网 络 图 的 结果 
来 看 , 单 倍 型 的 分 布 与 其 所 在 的 地 理 单元 之 间 不 存 
在 明显 的 对 应 关系 , 大 多 数 单 倍 型 此 处 于 较为 混 邓 
的 分 布 格局 中 。 

物种 遗传 多 样 性 的 高 低 与 其 进化 潜力 、 环 境 的 
适应 能 力 紧密 相关 , 较 高 的 遗传 多 样 性 对 于 物种 在 
环境 胁迫 的 情况 下 生存 繁衍 有 着 重要 的 意义 
( Hughes et al., 2008; Kearney et al., 2011; Wenzel 


小 车 蝗 的 线粒体 基因 具有 较 高 的 多 态 性 , 表明 该 种 
旦 忠 对 于 变化 的 环境 有 者 较 强 的 生存 能 力 , 这 可 能 
与 其 生物 学 特性 等 因素 有 关 。 李 春 选 等 (2003 ) 应 
用 等 位 酶 分 析 表 明 , 不 同 的 曝 虫 种 类 会 因 其 生物 学 
特性 及 对 栖息 环境 的 的 依赖 程度 , 体现 出 不 同 的 遗 
传 多 样 性 水 平 。 尽 管 黄 胜 小 车 蝗 一 般 锌 认为 是 分 布 
RNE, BAPER A Aaa RH 
的 选择 进行 定 回 扩散 , 生存 的 环境 复杂 多 样 ， 而 可 
能 正 是 生态 栖息 环境 的 多 样 性 , 造成 了 该 种 蝗虫 种 
群 内 丰富 的 遗传 多 态 性 。 田 外 , WENERA] 
地 理 维度 上 的 生态 适应 性 及 生活 史 的 多 样 性 , 也 可 
能 是 导致 种 群 遗 传 多 样 性 较 高 的 原因 。 

Tajima’ s D 检验 是 基于 种 内 多 人 态 性 的 一 种 中 性 
检验 方法 , 因 可 反映 出 物种 种 群 变 化 动态 的 历史 ， 
而 被 广泛 使 用 。 检 验 值 显著 大 于 0 时, 可 反映 出 种 
群 的 动态 平衡 ; 检验 值 显著 小 于 0 时 , 则 可 反映 出 
种 群 的 扩张 及 瓶 贷 效 应 (Tajima, 1989; Harpending 
et al., 1998)。 基 于 COI Fr EE BUR] ASSERERE BJ 
中 性 检验 结果 显示 , 总 群体 和 各 种 群 的 Tajima' s D 
AR EAS. 表明 序列 在 进化 上 遵循 中 性 模 
型 ,该 种 害虫 总 体 在 较 近 的 历史 上 未 经 历 群 体 扩张 
及 持续 增长 模式 , 保持 相对 稳定 的 群体 规模 。 

TRUE] AE EE BEES B] BS x RC Fst、 和 遗传 分 化 
系数 Gst 数值 均 较 低 , 基因 差异 程度 不 大 , 表明 各 
地 理 种 群 间 遗 传 分 化 程度 较 小 ( Rousset, 1997)。 从 
基因 流 Nm 人 研究 结 采 来 看 , 通常 认为 ，Nm «1 时， 
遗传 漂 变 成 为 种 群 间 遗 传 分 化 的 主要 因素 ,Nm >4 
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时 , 则 说 明 种 群 间 存 在 着 较 强 的 基因 流 , EAE 
化 作用 , 在 相当 程度 上 阻碍 了 种 群 间 的 遗传 分 化 
(Millar and Libby, 1991; Whitlock and Mccauley, 
1999)。 本 研究 中 , Xe E d ERES S) 3E EAT v 
Nm 29.01 >4, 表明 不 同 种 群 间 的 基因 交流 较 充 
分 。 分 子 变异 分 析 结 果 进 一 步 证 实 , 黄 胜 小 车 量 上 总 
群体 的 遗传 变异 来 主要 自 于 种 群 内 ,而 非 种 群 之 
间 , 反映 出 种 群 间 的 基因 交流 频繁 , 遗传 分 化 程度 
较 低 。 

黄 肥 小 车 蝗 种 群 间 相对 较 高 的 基因 流水 平 与 其 
生态 习性 密切 相关 。 郑 先 云 (2006 ) 人 研究 认为 , 虽然 
与 有 迁 飞 能 力 的 蝗 忠 相 比 较 , 黄 肥 小 车 蝗 种 群 间 基 
因 交 流水 平 较 低 , 但 却 高 于 所 人 研究 的 其 他 非 迁 飞 性 
蝗虫 。 相 对 于 稻 蝗 等 对 特定 生态 环境 依赖 性 强 、 分 
布 范 围 受 限 的 星 虫 种 类 ( 马 晋 等 , 2010) , RIDE 
尾 常会 因 环 境 的 变化 而 进行 迁移 。 呐 用 小 车 旦 的 曙 
晴 、 成 虫 都 善于 跳跃 ,成 虫 擅 飞 翔 。 据 河北 省 安 新 
县 植保 站 观察 , 黄 肥 小 车 蝗 3 龄 星 晴 无 植被 时 2 h 
迁移 10 m, 成 虫 会 因 作 物 生 长 情况 的 变化 而 发 生 较 
远 距离 的 迁移 ( 任 春光 , 2009)。 随 着 季节 的 变化 ， 
黄 肥 小 车 旦 从 草原 地 区 迁移 到 农 牧 交错 区 或 农田 进 
行 危 害 的 情况 也 较为 常见 。 可 能 正 是 这 种 个 体 或 者 
种 群 的 迁移 现象 , 增加 了 物种 内 的 基因 交流 ,弱化 
了 黄 胜 小 车 蝗 种 群 间 的 基因 分 化 。 男 一 方面 , 前 人 
经 研究 所 推测 的 星 总 科 进 化 呈现 出 的 简单 爆炸 式 辐 
HARA (Flook et al., 1995; Flook and Rowell, 1997; 
任 竹 梅 等 , 2003; A Bms, 2011). BRAR (MK 
JJ. KRE) 及 人 为 因素 的 影响 , 也 可 能 是 造成 黄 
肥 小 车 旦 地 理 种 群 间 发 生 基 因 交 流 的 原因 。 

本 研究 所 分 析 的 10 个 黄 胜 小 车 蝗 地 理 种 群 ， 
采集 于 7 个 省 份 , 有 着 较 大 的 地 理 距离 跨度 , HAB 
关 性 分 析 结 果 表 明 , 种 群 间 的 遗传 距离 与 地 理 距 离 
之 间 并 无 显著 相关 性 。 综 合 上 述 结果 可 知 ,， 黄 肥 小 
车 临 各 种 群 间 的 基因 交流 并 未 受到 地 理 距 离 的 影 
啊 , 地 理 隔离 并 非 是 导致 黄 盈 小 车 旦 种 群 间 遗传 分 
化 的 主导 因素 。 本 研究 工作 在 一 定 程 度 上 为 黄 肥 小 
车 晶 的 系统 演化 研究 提供 了 分 子 水 平 的 实验 证 据 ， 
为 该 种 害 忠 的 合理 有 效 控制 提供 分 子 生 物 学 方面 的 
基础 资料 。 尽 管 本 研究 有 意识 地 选取 进化 速率 较 快 
的 基因 序列 进行 分 析 , 但 所 选 的 序列 相对 较 短 ， 从 
人 研究 结果 来 看 , 变异 位 点 的 数量 较 少 , 还 不 足以 获 
得 全 面 的 遗传 学 数据 。 在 下 一 步 的 研究 中 , 还 需 增 
加 所 选 分 子 标记 测序 长 度 , 辅 以 其 他 分 子 标记 人 研究 
方法 , 增加 试验 采样 地 及 样本 数量 , 进一步 对 相关 


问题 进行 深入 分 析 ， 从 而 更 加 全 面 地 了 解 黄 胜 小 车 
星 各 分 布地 区 的 种 群 间 关 系 及 遗传 机 制 。 
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